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A "spin-flip transistor" is a lateral spin valve consisting of ferromagnetic source drain contacts to 
a thin-film normal-metal island with an electrically floating ferromagnetic base contact on top. We 
analyze the rfc-current-driven magnetization dynamics of spin-flip transistors in which the source- 
drain contacts are magnetized perpendicularly to the device plane by magnetoelectronic circuit 
theory and the macrospin Landau-Lifshitz-Gilbert equation. Spin flip scattering and spin pumping 
effects are taken into account. We find a steady-state rotation of the base magnetization at GHz 
frequencies that is tuneable by the source-drain bias. We discuss the advantages of the lateral struc- 
ture for high-frequency generation and actuation of nanomechanical systems over recently proposed 
nanopillar structures. 

I. INTRODUCTION 

Current induced magnetization excitation by spin-transfer torquei^ attracts considerable attention because of po- 
tential applications for magnetoelectronic devices. The prediction of current-induced magnetization reversal has been 
confirmed experimentally in multilayers structured into pillars of nanometer dimensionsA^£*& The devices typically 
consist of two ferromagnetic layers with a high (fixed layer) and a low coercivity (free layer), separated by a normal 
metal spacer. The applied current flows perpendicular to the interfaces. Often magnetic anisotropies force the mag- 
netizations into the plane of the magnetic layers. Recently a number of theoretical proposals pointed out interesting 
dynamics when the magnetization of one of the layers is oriented perpendicular to the interface planes Ji-tP i 

Fundamental studies of charge and spin transport have also been carried out in thin-film metallic conductors 
structured on top of a planar substrateii2iii*i2ii&i!ii£ The advantages compared to pillar structures are the flexible 
design and the relative ease to fabricate multi-terminal structures with additional functionalities such as the spin- 
torque transistor^ The easy accessibility to microscopic imaging of the structure and magnetization distribution 
should make the lateral structure especially suitable to study current-induced magnetization dynamics. Previous 
studies focused on the static (dc) charge transport properties, but investigations of the dynamics of laterally structured 
devices are underway^ii^ Recently, non-local magnetization switching in a lateral spin valve structure has been 
demonstrated^ In the present paper we investigate theoretically the dynamics of a lateral spin valve consisting of 
a normal metal film that is contacted by two magnetically hard ferromagnets. As sketched in Fig. a (nearly) 
circular and magnetically soft ferromagnetic film is assumed deposited on top of the normal metal to form a spin-flip 
transistor^ We concentrate on a configuration in which the magnetization direction of the source-drain contacts lies 
perpendicular to the plane of the magnetization of the third (free) layer. This can be realized either by making the 
contacts from a material that has a strong crystalline magnetic anisotropy forcing the magnetization out of the plane, 
such as Co/Pt multilayerSfSior by growing the source/drain ferromagnetic contacts into deeply etched groves to realize 
a suitable aspect ratio. In such a geometry, the magnetization of the free layer precesses around the demagnetizing 
field that a rises when the magnetization is forced out of the plane by the spin-transfer torque, as has been discussed 
in Refs. F^ l8l9l Therefore, as long as the out-of-planc magnetization of the free layer remains small, the free layer 
magnetization will always stay almost perpendicular to the source and drain magnetizations. In the present article we 
analyze in depth the coupled charge-spin-magnetization dynamics in such current-biased thin-film "magnetic fans" 
and point out the differences and advantages compared to the perpendicular pillar structures. A convenient and 
accurate tool to compute the dynamic properties of our device is the magnetoelectronic circuit theory for charge 
and spin transport^ coupled to the Landau-Lifshitz-Gilbert equation in the macrospin model. We include spin flip 
scattering in normal and ferromagnetic metals and the spin-pumping effect^S^ 

The article is organized as follows: In Section [D] we briefly review the Landau-Lifshitz-Gilbert equation including 
the current driven and spin-pumping torques that can be derived by circuit theory. In Section [fill the specific results 
for our "magnetic fan" are presented. The potential applications will be discussed in Scction llVl Section[V]is devoted 
to the conclusion. 
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II. FORMALISM 



We are interested in the magnetization dynamics of the soft ferromagnetic island (i.e., composed of permalloy) 
on top of the normal film as sketched in the Fig. ^ The Landau-Lifshitz-Gilbert (LLG) equation in the macro- 
spin model, in which the ferromagnetic order parameter is described by a single vector M with constant modulus 
M s , appears to describe experiments of current-driven magnetization dynamics well, 24 although some open questions 
remain Micromagnetic calculations of the perpendicular magnetization configuration in the pillar structure suggest 
a steady precession of the magnetization/^ The LLG equation for isolated ferromagnets has to be augmented by the 
magnetization torque L that is induced by the spin accumulation in proximity of the interface as well as the spin 
pumping: 

1 dm a dm 1 

— = -m x H ef f H m x — h — r L 1 

7 dt 11 7 dt VM S v ' 

where 7 is the gyromagnetic constant, m — M/Af s and H e ff is the magnetic field including demagnetizing, anisotropy 
or other external fields, ao is the Gilbert damping constant and V is the volume of the isolated bulk magnet. 

L = m x (iW+lP)) x m. 



where li P ^ and li b ' denote the pumped^ and bias-driveniS spin currents leaving the ferromagnet, respectively, and 
the vector products project out the components of the spin current normal to the magnetization direction. 

In magnetoelectronic circuit theory a given device or circuit is split into nodes and resistors. In each node a charge 
potential and spin accumulation is excited by a voltage or current bias over the entire device that is connected to 
reservoirs at thermal equilibrium or by spin pumping. The currents are proportional to the chemical potential and 
spin accumulation differences over the resistors that connect the island to the nodes. The Kirchhoff rules representing 
spin and charge conservation close the system of equations that govern the transport. In the following we assume 

that the ferromagnetic layer thickness is larger than the magnetic coherence length A c = ir/ k F — k F in terms of the 

majority and minority Fermi wave numbers that in transition metal ferromagnets is of the order of Angstroms. 
Let us consider a ferromagnet-normal metal (F\N) interface in which the ferromagnet is at a chemical potential 
and spin accumulation fifm (with magnetization direction m), whereas the normal metal is at and spin 
accumulation s. The charge current (in units of Ampere) and spin currents (in units of Joule), into the normal metal 



are2& 



Ic = [2g(Ho -»o)+ P9Vs ~ P9 m ■ s l ( 2 ) 

I<» = £lM^ ~ Vo) + Mf - (1 - Vr)m ■ s]m 

9 9 , / N 




FIG. 1: The model system consists of hard-magnetic source and drain contacts (Fl and F2) with antiparallel magnetizations 
perpendicular to the plane. On the top of the normal metal N, a soft ferromagnetic film (F3) is deposited with a slightly 
elliptical shape. The quantization direction, i.e., z-axis, is chosen parallel to the magnetization of the source and the drain. 
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where [j,q and are the chemical potentials in the ferromagnets and normal metal, respectively, g', g* are the 
dimensionless spin dependent conductances with polarization p — (g' — g^)/(g^ + g^) and total contact conductance 
g = g' + g*. In the Landauer-Biittiker formalism 

ff ™=M-5>?fo| a (4) 

nm 

where M is the total number of channels and r"™^ is the reflection coefficient from mode m to mode n for spin 
up(down) electrons. The spin transfer torque is governed by the complex spin-mixing conductance g*^, given by2& 

^=M-2»r(rr)*. ( 5 ) 

nm 

introduced in Eq. |JSJ in terms of its real and imaginary part as r\ r — 2Re<?^/g and r\i = 2lmg^/g. All conductance 
parameters can be computed from first principles as well as fitted to experiments. 
Slonczewski's spin transfer torque can then be written as 

m x x m = — ?7r[s — (s • m)m] + -^-rji(s x m). (6) 



The spin-pumping current is given by2£ 



I r ) = -zz9[rhinx—+r k —) . (7) 



8tt V dt dt 

We consider for simplicity the regime in which the spin-flip diffusion length l^j in the normal metal node is larger 
than the size of the normal metal region^ Charge and spin currents into the normal metal node are then conserved 
such thatSi 



i 

1^+1?])= I s /. (9) 



where we introduce a leakage current due to the spin-flip scattering I s / = g s fs/Ait and g s f — hvoosV 'n /t£j is the 
conductance due to spin flip scattering, where vdos is the (on-spin)density of state of the electrons in the normal 
metal, rj^ is the spin flip relaxation time and Vjv the volume of the normal metal node. 

The polarization of the source-drain contacts is supposed to be an effective one including the magnetically active 
region of the bulk ferromagnet with thickness governed by the spin-flip diffusion length in the ferromagnet. For the 
free magnetic layer F3, the perpendicular component of the spin current is absorbed to generate the spin transfer 
torque. The collinear current has to fulfill the boundary conditions in terms of the chemical potential /if = /if — /i| 
governed by the diffusion equation 

d 2 vf(z) _ fJ.f(z) ^ 



dz 2 

where lf d is the spin flip diffusion length in the ferromagnet 



III. SPIN TRANSFER TORQUE AND MAGNETIC FAN EFFECT 

In this Section, we solve the Landau-Lifshitz-Gilbert equation including expressions for the spin-transfer torque on 
the free layer according to the circuit theory sketched above. 

A. Currents and spin torque 

In metallic structures the imaginary part of the mixing conductance is usually very small and may be disregarded, 
i.e., rji ~ 0. The source and drain contacts FX\N and F2\N are taken to be identical: g\ = g2 = g, Pi = P2 = P and 
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Vri = Vr2 = i)r- For .F3|iV we take rj r 3 = 773. In our device, the directions of the magnetization of the fixed magnetic 
leads are ni! = (0,0, 1) and m 2 = (0,0, —1). For the free layer we allow the magnetization m 3 = (m x ,m y ,m z ) to be 
arbitrary. We assume that F3 is a floating contact in which the the chemical potential [j,q 3 adjusts itself such that 
the net charge current through the interface F3\N vanishes: 

4 3) = |^[2(Mo F3 - tf) +P3Mf 3 -Pss ■ m 3 ] = 0. (11) 



Applying a bias current Jo on the two ferromagnetic leads, JT and F2, the conservation of charge current in the 
normal 1 
dictates 



normal metal then gives iP = -J< 2) = Jo. At the F3\N interface, the continuity of the longitudinal spin current 



where <J-\{<j{) is the bulk conductivities of spin up (down) electrons in the ferromagnet and A the area of the interface. 
Choosing the origin of the z axis is at the F3\N interface and assuming F3 to be of thickness d, 

dm \ ( dn i 



"^J--"^J--°- (13) 

With both boundary conditions, the diffusion equation can be solved for the spin accumulation in FZ 

C 3 cosh(^)s- m 3 

(*) = T —1 ( 14 ) 



C 3 + a tanh( A- ) cosh( -4- ) 

sd J sd 

where £ 3 = 53(1 — pi)/4 characterizes the contact F3|N and a = hAa^a \ / '(e 2 lf d (a^ +crj.)) describes the bulk conduction 
properties of the free layer with arbitrary m 3 . The limit d -C lf d corresponds to negligibly small spin-flip, which implies 
tanh (d/lf d ) ~ 0. Near the interface, the spin accumulation in F3 then reduces to 

/if 3 = s ■ m 3 . (15) 

(3) 

In this limit, l s ■ m 3 = the collinear component of the spin current vanishes. 

By solving the linear equations generated by Eqs. I|8I9|I . we obtain the spin accumulation s in the normal metal 
node, 

s = C ■ [8tt4 p) + W b ] (16) 

where the elements of the symmetric matrix C are given in Appendix |X|and W b = (0, 0, 2phl /e) is a bias-vector. Eq. 
dim contains contribution due to bias current and spin pumping effect. The spin accumulation in the ferromagnet 
Eq. i|14|) should be substituted in Eq. (|16|) to give the spin accumulation in the normal metal, from which the spin 
transfer torque can be determined according to Eq. ©. For an ultrathin film, the spin transfer torque, including 
pumping effect and spin accumulation in the ferromagnet, reads, 



L = 11393 
8tt 

with the elements of II listed in Appendix. 



II • [8vrlW + W 6 ] , (17) 



B. Dynamics of the free layer 

After the bias current is switched on, a spin accumulation builds up in the normal metal. At the beginning, the 
spin-transfer torque exerted on the magnetization of the free layer (F3) causes a precession out of the plane, hence 
generating a demagnetizing field that is oriented perpendicular to the film plane. Subsequently the magnetization 
precesses around and as long as the current Jo continues, the rotation persists. In order to determine the dynamics 
of the magnetization, we apply the spin torque term L [Eq. lfT7| ] to the Landau-Lifshitz-Gilbert (LLG) equation 
Crystalline anisotropics in F3 may be disregarded for soft ferromagnets such as permalloy. The effective field in the 
LLG equation then reduces to 

H A = -^ M s (N x m x , Nytriy, N z m z ) , (18) 
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where N x , N y and N z are the demagnetizing factors determined by the shape of the filmiSi The anisotropy field 
keeps the magnetization in the plane when the torque is zero. The spin torque generated by the current bias forces 
the magnetization out of plane, hence triggering the nearly in-plane rotation of the magnetization. Substituting the 
spin-torque term Eq. H17fl into Eq. we obtain for the following LLG equation, 

1 dm 1 , dm „ , T „ . . 

-— = -m x H A + - (ao + a ') m x — + H flt (J ) (19) 
7 at 7 at 



Here the last vector 



H s t(/ ) = ^^stj^y(-m x m z ,-m y m z ,l - ml) . (20) 



is the effective field induced by the spin-transfer torque that depends on the position of the magnetization and the 
device parameter 

A = pnmGi m) 

St Gtg 3 + 2(p^-l + V )gg i (l-ml)' { ' 

where Gi's are introduced in Appendix A. According to Eq. 1|21[) . we can accurately engineer the device performance by 
tuning the conductances and polarizations. Compared with the original LLG equation, a new dimensionless parameter 
entering the calculation 

<_^. 7 S(Reo u ) 2 - , , 

2wVM s v ; 

reflects the tensor character of the pumping-induced additional Gilbert damping^ Choosing contact F3|iV to be 
metallic and the others to be tunneling barriers, the condition 173 3> g,g s f can be realized. In that limit ^ct' reduces 
to 

Regl 1 , (23) 



4ttVM l 

which agrees with the enhanced Gilbert damping derived in Ref. |23- In the following, we take a = ao + a' to be the 
enhanced Gilbert damping constant. 



1. Vanishing in-plane anisotropy 

Here we rewrite the free layer magnetization in two polar angles <j) (in-plane) and 9 (out-of plane) such that 
m = (cos 6 cos 0, cos 9 sin 0, sin 6) and assuming a small z-component, i.e., m z = sin# w 9 and cos0 w 1. When the 
free layer is a round flat disk with demagnetizing factors N x = N y s» and N z « 1, the Eqs. (|19|l reduce to: 

d<t> d9 7\ 7" f\ 

d- t = a i + ^> ^ 

introducing = hA s tIo / (2eM s V) . Eq. (|24(l separates the motion for the in and out-of-plane angles. We consider 

the dynamics of a current that is abruptly switched on to a constant value Iq at t = 0, assuming that 9{t — 0) = 0, 
i.e., a magnetization that initially lies in the plane. The motion of 9 for t > is then given by 



dO 



dt 1 + a 2 

where we introduced the response time 



a -tj 



0J^~ tlT . (25) 



< 1 + ° 2 » ,26) 



aptoiMsNz 
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In-planeRotation 




FIG. 2: The in-plane rotation (in the unit of giga hertz) versus time (in nano seconds). The solid line: including spin pumping 
effect. The dash line: without spin spumping effect. 



and the saturation in-plane rotation frequency 



(27) 



Similarly, the in-plane rotation is governed by 



cf>(t) = —uijt + 



jafj, M s N z 
d 4 = -U i + J^e-" 



dt 



l + a 2 



(28) 



Taking the parameters from Ref. Il2l viz. a volume of normal metal V n — 400 2 x 30 nm 3 , spin flip time in the 
normal metal of r s / = 62 ps, density of states v^os == 2.4 x 10 28 eV _ m -3 , we find e 2 g s f fh = 0.3 

Let us take the thickness of the free layer d = 5 nm. The saturation magnetization of permalloy is M s = 8 x 
10 5 A m . The relative mixing conductance is chosen 773 ~ r\ r ~ 1 and the bulk value of the Gilbert damping constant 
for Py is typically ao = 0.006. 22 A metallic interface conductance (for F3|iV) is typically 1.3JD m 1 ^, whereas the 
source/drain contacts are tunneling barriers with resistance hj (e 2 g) — 20 kflU The calculated enhancement of the 
Gilbert damping constant is then a' = 0.004 and the response time r = 0.52 ns. The motion of the magnetization of 
the free layer is depicted by Fig. [21 and Fig. for a bias current density of J — 10 7 A cm -2 with the cross section at 
the electronic transport direction 400 x 30 nm 2 iS 

The spin pumping effect through the enhanced Gilbert damping constant reduces the saturation frequency from 
2.0 to 1.2 GHz , but also the response time to reach the saturation value from 0.87 to 0.52 ns. Notice that the 
frequency is directly proportional to Jo and thus in the absence of any in-plane anisotropy the frequency can be tuned 
continuously to zero by decreasing the bias current. The out-of-planc motion is very slow compared to the in-plane 
one: it decreases from 12 MHz to around when the in-plane rotation approaches the saturation frequency. As shown 
in Fig. within a long period the small angle approximation still holds. A larger ratio of gs/g also gives higher 
frequencies. Decreasing the diameter, and thus also the volume, of the free layer gives a smaller demagnetizing factor 
N z , which causes larger a response time r according to Eq. (|26|l and increases the saturation value of the in-plane 
rotation frequency w^. 



2. In-plane anisotropy 



In reality, there are always residual anisotropics or pinning centers. Shape anisotropies can be introduced inten- 
tionally by fabrication of elliptic F3 discs. We consider the situation in which the free layer is slightly pinned in the 
plane by an anisotropy field that corresponds to an elliptic (pancake) shape of the ferromagnet. At equilibrium, the 



7 




FIG. 3: The out-of-plane motion(in the unit of mega hertz) versus time (in nano seconds). The solid line: including spin 
pumping effect. The dash line: without spin pumping effect. 



Out-of-plano^ngle 




FIG. 4: The out-plane angle 9 (in degree) versus time (in nano seconds). The solid line: including spin pumping effect. The 
dash line: without spin spumping effect. 



F3 magnetization is then aligned along the easy, let us say, x— axis. The in-plane rotation can be sustained only when 
the spin transfer torque overcomes the effective field generated by the shape anisotropy, hence a critical current I c 
for the steady precession is expected. For an ellipse with long axis of 200 nm, thickness 5 nm and aspect ratio 0.9, 
the two demagnetizing factors are calculated to be N y = 0.0224 and N x = 0.0191. With a Gilbert damping constant 
a = 0.01, the numerical simulation gives I c = 4.585 mA corresponding to a current density J c = 3.8 x 10 7 A cm~ 2 
(the cross section is 400 x 30 nm 2 )^ 

These critical current densities are of the same order of magnitude as those used to excite the magnetization in 
spin- valve pillars. So even a relatively small anisotropy can cause a significant critical current. In order to operate 
the magnetic fan at small current densities, the magnetic island should be fabricated as round as possible. The 
magnetization responds to a current step function below the critical value by damped in-plane and out-of-plane 
oscillations and comes to rest at a new in-plane equilibrium angle <j) e with zero out-plane component (c/. Figs. [5] 
and EJ) . At the steady state, the spin-transfer torque is balanced by the torque generated by the in-plane anisotropy, 
i.e. the angle <p e is determined by sin(20 e ) = 2!F(Iq )/{hqM s {N v — N x )). With given bias current, smaller \N y — N x \ 
correspond to larger in-plane angles \<p e \- According to the theory of differential equations^ the frequency for the 
damped magnetization oscillation can be found by diagonalizing the LLG equation at the "equilibrium point" given 
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X-component of Magnetization 




FIG. 5: Below critical current, the x-component of magnetization versus time (in nano seconds). The bias current is 4.5 mA. 
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FIG. 6: Below critical current, the ^-component of magnetization versus time (in nano seconds). The bias current is 4.5 mA. 



by 



this leads to 



l^yJ(2N z -N x -N y )y/v(lo)+V(I ) , 



(29) 



where 



»«■>-<*.-*.>■- ^ 



(30) 



Equation l|29|l teaches us that below the critical current, decreasing the current increases the rotation frequency. 
Changing the damping constant does not change ui^ for a given current but only changes the response time to reach 
the new equilibrium. 

As shown by Fig. to Fig.[5|the magnetization above the critical current saturates into a steady precessional state 
accompanied by an oscillation of the z-component (nutation) . In this situation, <p e is no longer a constant of motion. 



Instead the new steady state is given by m 2 



= and m z — !F{Iq) / '(afJ,oM a N z ). Diagonalizing the LLG around 



this point we derive the in-plane rotation frequency 



y/(N z -N x )(N z -N y ) 



n 



(31) 



In the limit of vanishing in-plane anisotropy, i.e., N x = and N y = 0, we recover the previous result. As shown by 
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X-component of Magnetization 




FIG. 7: Above the critical current, the ^-component of magnetization versus time (in nano seconds). The bias current is 
4.6 mA. The frequency is about 3.6 GHz. 
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FIG. 8: Above critical current, the ^-component of magnetization versus time (in nano seconds). The bias current is 4.6 mA. 



Trajectory of Magnetization 




FIG. 9: Above critical current, the trajectory of magnetization within 5 nano seconds. The bias current is 4.6 mA. This picture 
clearly shows the steady precession of the magnetization. 
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FIG. 10: The critical current I c versus damping constant a. This figure shows saturation of 7 C above a critical a. 



Fig. 1101 the dependence of the critical current on the damping constant is different from the simple proportionality 
predicted for pillar structures^ Specifically we observe saturation of the critical current above a critical damping. 

In the anisotropic case the extra power necessary for maintaining the motion generates microwaves which may 
be attractive for some applications. 



IV. APPLICATIONS 



Our "magnetic fan" has the advantage that the magnetization dynamics is not hidden within the structure as in 
the pillars, but is open to either studies of the dynamics by fast microscopy, or to the utilization of the dipolar field 
from the soft magnetic island. We envisage applications as magnetic actuators for nanomechanical cantilevers and 
nanoscale motors, as nanoscale mixers of biological or biomedical suspensions containing magnetic nanoparticles, or 
as magnetic resonance detectors, again possibly useful for biomedical applications. 



A. Actuators 



The rotating magnetization of the "magnetic fan" generates a periodic dipolar field which can be applied to actuate 
a nanomechanical cantilever with a (hard) ferromagnetic tip. Assuming for simplicity that the magnet F3 and the 
cantilever are at a sufficiently large distance the force on the cantilever magnet is given by 

F = KV(M C • Hi) , (32) 

where M c is the saturation magnetization and V c is the volume of the cantilever magnet and the field Hd generated 
by a magnetic dipolar at the position r can be written as 

H,^ g±^ . (33) 

Assume a cantilever on top of the magnetic fan at a distance of 125 nm (along ^-direction),-! with beam plane parallel 
to the plane of the Py film F3 and magnetization along the x-axis. The saturation value of cantilever magnetization 
is taken as 1.27 x 10 6 A m _1 . Assuming a lateral size of the cantilever magnet of 150 x 150 nm 2 with thickness 
50 nm, the force is estimated to be 

F = 1.1 X 10~ 8 cos(cV) N (34) 



where oj^ is the rotation frequency of the "magnetic fan" . To efficiently generate the mechanical modes of the 
cantilever, the cantilever magnet should be hard enough. 
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Fixing other parameters, the force scales like 1/r 4 with respect to distance r. When the two ferromagnets are closer 
to each other the distribution of the magnetizations increases the force over the value estimated above. We see that 
in the dipole-approximation, the force is already quite significant and it will be significantly larger when the the full 
magnetostatic energy is computed. 

Generally, the torque on the cantilever may generate both flexural and torsional motion on the cantilever. The 
torsional motion coupled to the magnetization dynamics has been investigated for such a system 3 ^ and the nanome- 
chanical magnetization reversal based on the torsional modes has been proposed* 3 ^ The coupling of a cantilever to 
the oscillating dipolar field will be discussed elsewhere. 

B. Mixers 

The dipolar field produced by our device can also be used to function as mechanical mixer for suspensions of 
magnetic particles. To this end we should scale down the frequency of the rotating magnetization either by decreasing 
the bias current or re-engineering the parameters of the device, e.g., increasing the thickness of the Py film. Low 
saturation magnetization is detrimental in this case, since that would also reduce the usable stray fields. By these 
ways, one hopefully can access the kilo hertz frequency region, which is important for the hydrodynamic motion in 
ferrofluidsi^ 



C. Detectors 

An external field influences the frequency of the rotation of the magnetization. Response to the change of the 
frequencies is the rebuilding of the spin accumulation in the normal metal hence altering the source-drain resistance 
Rsd- Due to the relation 

^ - Mo 2 = RsdIo , (35) 

this deviation is reflected on the source-drain voltage-current curve. This feature can be implemented as a sensor for 
biomedical applications in order to detect the presence of magnetic beads, which are used as labels in biosensors: 35 
Furthermore, the ability to change the frequency of the "magnetic fan" should allow to measure locally the frequency 
dependence of the magnetic susceptibility, which offers an alternative pathway to using magnetic nanoparticles for 
biosensing applicationsi 36 i 37 

V. CONCLUSION 

We studied the magnetization dynamics of a magnetic transistor, i.e., a lateral spin valve structure with 
perpendicular-to-plane magnetizations and an in-plane free layer attached to the normal metal that is excited by 
an external current bias. By circuit theory and the Landau-Lifshitz-Gilbert equation, analytic results were obtained 
for the spin-transfer torque and the dynamics of the magnetization in the limit of small out-of-plane angle 9. Spin 
flip and spin-pumping effects were also investigated analytically and an anisotropic enhanced Gilbert damping was 
derived for the free layer magnetization. Without an externally applied magnetic field, a continuous rotation of the 
magnetization of the free layer at GHz frequencies can be achieved. In the lateral geometry, the free layer is no longer 
buried or penetrated by a dissipating charge current, thus becomes accessible for more applications. Our methods 
handle the microscopic details on crucial issues like spin-torque transfer efficiency, spin-flip scattering and the response 
time, hence offering accurate design and control. The rotation can be observed, e.g., by magneto-optic methods. This 
new device has potential applications as high frequency generator, actuator of nanomechanical systems, biosensors, 
and other high-speed magnetoelectronic devices. 
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APPENDIX A: SPIN ACCUMULATION IN NORMAL METAL NODE 

Here we summarize a number of complex angle dependent coefficients. The elements of the symmetric matrix C in 
Eq.JTJJ read 

n Gt(Gi - Gimp - 2g{p 2 -1 + V )(g t - G 4 m 2 y ) 

On — (Al) 

„ ri GiG±rn x m, y GtGim x m z 

C12 =C 2 i = -, and Ci 3 = C 3 i = (A2) 

n _ G2{Gt - Gimp - GtGjm 2 

22 — q (AO) 

C 23 =C 32 = ^» and C 33 = Gt{Gl + Q Qiml) (A4) 

introducing 

Si = (1 - P3XI - A 3 ).93 + 2r7.g + 2g sf (A5) 

G2 =m93 + 2(l-p 2 )g + 2g sf (A6) 

& = (1 - Pl)(l - A 3 )g 3 + 2(1 - p 2 )g + 2g sf (A7) 

e?4 = r?393-(l-P3)(l-A 3 )ff 3 (A8) 

St = %.93 + 2?7<7 + 2 3s/ (A9) 

Q = Gt[GtG 3 + 2{ P 2 - 1 + v)gGi(i - m 2 )] (Aio) 

C 3 + c tanh(d/Z; d ) 

in the limit of negligible spin flip in F, i.e., d -C Z^, then A3 rj 1. The elements of the matrix in Ea. (|17|l are given by 



g t g 3 (l - mp + 2Q^p 2 - 1 + V )gmz 
Q 



n n xJ "» y (A12) 



„ „ -GiG2m x m y -GtGim x m z 

II12 =II 21 = 2-, and II13 = (A13) 

_ StS 3 (l-ml) + 2g 4 (p 2 -l + V )gml 

n 22 = v - (Ai4) 

„ -GtGim y m z -g t g 3 m x m z 

n 2 3 = g- 2 , and n 3i = (A15) 

„ -GtGzm,ym z GtGi(l-m 2 z ) 

n 32 = g- 2 , and n 33 = (A16) 

(A17) 
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